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HIGHLIGHTS 


•  An  adaptive  curve  transformation  is  used  to  quantify  the  battery  capacity  fade. 

•  Only  one  “time-based”  parameter  is  needed  in  the  model  to  enable  SOH  estimation. 

•  Only  the  initial  charging  data  is  needed  to  train  the  estimation  model. 

•  The  SOH  estimation  accuracy  is  not  affected  by  duty  cycles  or  aging  history. 
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An  effective  method  to  estimate  the  state  of  health  (SOH)  of  lithium  ion  batteries  is  illustrated  in  this 
work.  This  method  uses  an  adaptive  transformation  of  charging  curves  at  different  stages  of  life  to 
quantify  the  extent  of  capacity  fade  and  derive  a  time-based  parameter  to  enable  an  accurate  SOH 
estimation.  This  approach  is  easy  for  practical  implementation  and  universal  to  chemistry  or  cell  ge¬ 
ometry,  with  minimal  demand  of  learning.  With  a  typical  constant  current— constant  voltage  (CC— CV) 
charging  method  for  a  lithium  ion  battery,  this  approach  uses  an  equivalent  circuit  model  to  characterize 
the  CC  portion  of  the  charging  curve  and  derive  a  transformation  function  and  a  time-based  parameter  to 
estimate  SOH  at  any  stage  of  life  via  a  nonlinear  least  squares  method  to  identify  model  parameters.  The 
SOH  estimation  errors  (discrepancy  between  estimated  and  experimental  values,  denoted  as  ASOH)  are 
under  2%  before  the  end  of  life  in  cases  shown  at  25  °C  and  60  °C  and  a  range  of  typical  discharging  rates 
up  to  3C.  With  different  sizes  and  chemistries,  the  ASOHs  are  all  less  than  3%. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Lithium  ion  batteries  (LIBs)  are  widely  used  in  portable  elec¬ 
tronic  products  and  new  energy  technologies  such  as  electric  ve¬ 
hicles  and  electrical  energy  storage  systems.  They  have  some 
incomparable  advantages  compared  to  other  rechargeable  batte¬ 
ries,  including  higher  energy  density,  less  maintenance  require¬ 
ment,  higher  nominal  voltage,  lower  self-discharge  rate;  thus,  LIB  is 
regarded  as  one  of  the  most  promising  energy  solutions.  However, 
its  performance  degrades  with  time  because  of  impacts  from 
operating  environment  and  aging  history.  The  battery  failure  may 
result  in  serious,  sometimes  catastrophic,  consequences.  Therefore, 
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it  is  desirable  to  estimate  the  state  of  health  (SOH)  of  the  battery  to 
predict  how  soon  it  will  fail  or  void  the  guarantee  of  satisfactory 
performance. 

Extensive  studies  have  been  carried  out  for  the  SOH  estimation 
in  recent  years.  Some  consider  SOH  estimation  as  a  black-box 
problem,  using  artificial  intelligence  (AI)  and  machine  learning 
(ML)  to  estimate  SOH  for  various  operating  conditions  without  full 
understanding  of  the  aging  mechanism.  AI  and  ML  algorithms 
including  support  vector  machine  (SVM),  probabilistic  neural 
network  (PNN),  recurrent  neural  network  (RNN),  and  fuzzy  iden¬ 
tification  have  been  employed  in  various  attempts  for  SOH  esti¬ 
mation  [1—5].  However,  a  sufficient  amount  of  test  data  must  be 
obtained  to  train  the  estimation  model,  and  proper  estimation 
methods  and  algorithms  are  crucial  for  implementation.  There  are 
other  attempts  to  develop  numerical  models  to  describe  battery 
aging  behavior  in  mathematical  formulae,  equivalent  circuit 
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diagrams,  or  electrochemical  equations  [6-9].  In  these  models, 
specific  parameter  identification  techniques  are  used  to  correlate 
electrochemical  properties  in  the  LIB  with  degradation,  and  the 
SOH  is  reflected  by  one  parameter  or  combinatory  variants  deter¬ 
mined  using  parameter  identification  techniques.  In  Refs.  [6,7], 
simple  equivalent  circuit  models  (ECM)  were  proposed  to  support 
parameter  identification  methodology.  The  internal  resistance  and 
open  circuit  voltage  (OCV)  were  used  in  the  parameterization.  In 
Refs.  [8,9],  electrochemical-based  single  particle  model  (SPM)  were 
proposed,  and  a  nonlinear  least  squares  method  was  used  for 
parameterization.  These  approaches  mentioned  above  mostly 
derive  qualitative  SOH  estimation  based  on  state  of  charge  (SOC), 
impedance,  or  OCV,  which  is  not  as  user-friendly  as  expected  in 
actual  operations. 

In  this  paper,  an  effective  SOH  estimation  method  for  single  cells  of 
LIB  is  proposed.  The  same  approach  can  be  applied  to  multi-cell  sys¬ 
tems,  although  it  is  not  demonstrated  here.  This  approach  can  facili¬ 
tate  easier  implementation,  less  demand  for  prior  learning,  and  better 
generalization  to  accommodate  variations  in  cell  design  and  chem¬ 
istry.  Some  important  aspects  of  this  approach  are  denoted  here: 

1.  The  voltage-time  data  in  the  constant-current  (CC)  step  of  the 
charging  profile  are  adopted  in  this  approach  due  to  the 
following  considerations:  (1)  in  practical  use,  the  charging 
process  is  executed  in  a  more  controlled  manner  in  contrast  to 
discharging  process  under  a  duty  cycle;  (2)  the  voltage-time 
variables  are  easy  to  measure  than  the  others  such  as  resistance, 
OCV,  SOC  et  al.;  and,  (3)  the  constant- voltage  (CV)  step  in  the 
charging  profile  is  not  considered  because  it’s  more  time 
consuming  and  inconvenient  in  data  collection  and  analysis. 

2.  To  deduce  SOH  estimation  mathematical  model,  the  voltage 
variations  with  time  in  the  CC  step  of  the  charging  profile  are 
characterized  as  a  function  of  aging  duration.  A  typical  ECM  used 
by  many  others  is  introduced  to  derive  aging  parameters  for  the 
SOH  estimation.  Among  all  necessary  model  parameters,  only  a 
“time-based”  one  is  used  in  the  model  formula  to  enable  SOH 
estimation;  avoiding  the  use  of  OCV  or  SOC  as  variables  in  the 
estimation  process,  due  to  the  difficulty  in  deriving  these  vari¬ 
ables  in  practical  use.  This  unique  aspect  offers  very  attractive 
and  favorable  operability  in  actual  implementation. 

3.  Only  the  initial  charging  curve  as  a  benchmark  of  a  fresh  battery 
to  represent  the  chemistry  in  the  cell,  corresponding  to  the  100% 
SOH,  is  needed  in  this  model  and  the  initial  model  parameter¬ 
ization  is  learned  with  regression  algorithm  from  the  voltage¬ 
time  offline  data  in  the  CC  step.  With  the  aging  of  the  battery, 
the  SOH  can  be  quantified  by  the  unique  time-based  parameter 
in  the  estimation  model  formula,  and  the  parameters  identifi¬ 
cation  is  realized  by  a  nonlinear  least  squares  method. 

4.  As  long  as  the  operating  limits  of  voltage,  current,  and  tem¬ 
perature  are  not  exceeded,  the  voltage  profile  with  the  same 
charging  current  can  reflect  the  aging  of  the  battery  and  the 
variation  of  such  voltage  profile  with  aging  can  be  quantified 
through  a  continuous  update  of  the  time-based  parameter.  In 
other  words,  the  accuracy  of  the  estimation  model  is  not 
affected  by  duty  cycles  or  complicated  aging  history;  thus,  the 
applicability  of  this  SOH  estimation  method  is  quite  universal 
and  very  easy  to  implement. 


2.  Mathematical  model  for  SOH  estimation 

2.1.  Model  description 


Cycle  number 


Fig.  1.  (a)  Charging  curves  of  a  LiCo02/LiNixMnyCoz02-based  LIB  at  100%,  95%,  90%,  85%, 
and  80%  SOHs.  Typical  CC-CV  charging  method  is  used  and  only  the  CC  profiles  at  C/2 
are  shown,  (b)  The  corresponding  SOH  curve  as  a  function  of  cycle  number. 


profiles  were  recorded  using  a  typical  CC-CV  charging  method,  but 
only  the  CC  profiles  are  shown.  Fig.  1(b)  exhibits  the  SOH  changes 
with  cycle  number.  Intuitively,  these  charging  curves  are  similar  to 
one  another,  progressively  shifting  to  a  shorter  time  span  and 
higher  voltage  range  with  cycle  number.  Such  a  trend  inspires  this 
approach  of  SOH  estimation. 

These  CC  charging  profiles  can  be  represented  properly  with  an 
ECM  [6,7,10]  (Fig.  2).  The  ECM  consists  of  an  ohmic  resistor  (Rq),  a 
parallel  RC  network  (Rp  and  Cp),  and  a  voltage  source  ( Uoc )•  The  Rq 
represents  a  physical  resistive  effect  in  the  dynamics  of  the  battery 
kinetic  behavior,  which  is  responsible  for  the  instantaneous  IR 


Rq  Rp 

i — WV — i — WV — r^+ 


Uoc(SOC) 


u 


Fig.  1(a)  shows  the  charging  profiles  of  a  LiCo02/LiNixMnyCoz02- 
based  LIB  at  100%,  95%,  90%,  85%,  and  80%  SOH,  respectively.  These 


Fig.  2.  A  simple  ECM  of  LIB. 
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voltage  drop  in  the  charging  or  discharging  process  upon  polari¬ 
zation.  The  parallel  ftp-Cp  network  describes  the  nonlinear,  faradic 
charge-transfer  impedance  of  a  battery.  Uoc{SOC )  is  the  equilibrium 
OCV  of  the  battery,  increasing  nonlinearly  and  monotonously  with 
SOC.  U  is  the  battery  terminal  voltage  and  /  the  current  under  po¬ 
larization.  Two  assumptions  associated  with  ECM  are  given  below 
to  simplify  the  illustration  in  the  discussion: 


practice.  This  aspect  solved  the  dilemma  often  experienced  in  BMS 
where  the  Uoc{SOC )  is  difficult  to  measure  or  estimate.  However,  in 
Eq.  (6)  another  dilemma  in  dealing  with  SOC  remains.  SOC  is 
difficult  to  handle  because  the  parameter  cannot  be  determined 
directly  or  estimated  easily  with  certainty.  To  resolve  this  issue, 
Eq.  (2)  shall  be  used  to  represent  SOC  in  Eq.  (6).  Eq.  (6)  can  now  be 
rewritten  as 


1.  Rq,  ftp,  and  Cp  are  invariant  in  the  operating  range  (e.g.  3. 0-4.2  V), 
and  they  change  with  aging  of  the  battery. 

2.  Uoc(SOC)  is  only  associated  with  SOC  within  certain  error 
tolerance  [11  ].  The  Uoc(SOC )  curve  is  also  considered  invariant  in 
the  temperature  range  and  aging  process  during  operation. 

During  CC  charging  the  battery  voltage  can  be  expressed  by  a 
nonlinear  function  with  variables  SOC  and  time  t  as 

U{SOC,t)  =  Uoc(SOC)  -  Up-exp-5  +  /-Rp(l  -  exp-flj  +  IRa 
=  UodSOC)  -  {Up  +  l-Rp)exp---  +  /•  (Rp  +  RQ) 


U'(t')  =  U(t)  +  {Up  +/  Rp)exp-^t  -  (u'p  +  I-R'p)exp-7  t' 

+  l(R'p  +  R'p)  -/(Rp  +  Rfl)  (7) 


where  t!  corresponds  to  t  at  the  same  SOC,  but  at  a  different  age. 

If  the  U(t)  is  obtained  from  a  fresh  battery  at  SOH  =  100%,  Cs  ~  0, 
the  SOH  estimation  can  be  deduced  from  Eqs.  (7)  and  (8)  through  a 
further  simplification  as  follows:  by  making  k  =  CJC u,  At  =  C^/J, 
ci  —  Up  +  /•  ftp,  b  —  —  1/t,  c  =  Up  +  /  ftp,  d  =  —  1  /t  ,  e  =  /(ftp  +  Rq  —  Rp  —  Rq) 
( 1 )  a  universal  mathematical  expression  is  derived  as 


Where  z  =  Rp  •  Cp,  Up  the  polarization  voltage  when  t  =  0,  and  all 
variables  positive. 

In  Eq.  (1),  the  relation  between  SOC  and  t  can  be  further  estab¬ 
lished  by  introducing  the  recharged  capacity  Cr,  stored  capacity  Cs, 
and  usable  capacity  Cu.  The  Cr  represents  the  capacity  obtained 
from  the  CC  charging  step  started  from  t  =  0.  The  Cs  represents  the 
remaining  capacity  after  the  last  discharging,  and  the  Cu  represents 
the  actual  or  usable  capacity  as  the  results  of  aging.  All  of  above  can 
be  described  by  Eq.  (2) 


SOC  = 


Cr  +  Cs  _  /-t  +  Cs 

Cu  Cu 


Accordingly,  Eq.  (1)  can  be  rewritten  as 


(2) 


U{SOC,  t )  =  Uoc{SOC)  -  {Up  +I-Rp)exp-%A±^±%r51  +  I-{RP  +  Rq ) 

(3) 

The  battery  voltage  as  a  function  of  SOC  can  be  expressed  as 


U(SOC)  =  U0dSOC)-(Up  +  IRp)exp-&socexp5  +  I(Rp  +  RQ) 

(4) 


with  the  cycle  number  increasing,  Cu  will  degenerate  gradually.  At  a 
certain  state  of  aging,  Cu  is  denoted  as  C'u.  Accordingly,  the  pa¬ 
rameters  U,  Up,  Rp,  Rq,  Cs,  and  z  in  Eq.  (4)  are  deformed  as  U',  Up,  ft'p, 
R'Q,  C',  and  z.  On  the  premise  of  assumptions  1  and  2,  Eq.  (4)  can  be 
rewritten  to  describe  the  change  of  the  battery  voltage  as 


LZ'(t')  =  U(k(t'  +  At))  +  a-expbk(t+At)  -  c  expd  t’  +  e  (9) 

The  parameters  {k,  At,  a,  b,  c,  d,  and  e )  in  Eq.  (9)  can  be  identified 
by  using  a  numerical  calculation  method  to  yield  an  estimation  of  k. 
The  SOH  can  then  be  represented  as  1  Ik. 

2.2.  Parameter  identification  procedure 

As  a  LIB  ages,  its  capacity  fades,  as  indicated  by  the  gradually 
shortened  span  in  the  charging  profile  under  the  same  charging 
condition  and  the  rising  cell  voltage  at  the  same  capacity  as  shown 
in  Fig.  1(a).  This  gradual  change  in  the  charging  profile  as  a  function 
of  battery  aging  makes  it  possible  to  use  the  change  in  voltage¬ 
time  curves  to  estimate  SOH  and  battery’s  remaining  service  life. 

To  perform  parameter  identification,  we  shall  start  with  the  first 
charge/discharge  cycle  of  a  fresh  LIB  cell  as  a  baseline.  The  CC 
charging  profile  is  denoted  as  ( U ,  t).  By  using  a  nonlinear  least 
squares  method,  a  single-variable  function  U{t)  is  obtained  to 
represent  this  CC  charging  profile  and  used  it  as  a  benchmark  for 
100%  SOH.  The  parameters  a  and  b  should  be  identified  first.  In  the 
parameterization  process,  several  constraints,  including  k  -»  1, 
A t  -►  0  and  e  -►  0,  must  be  met.  Once  that  a  ~  c  and  b  ~  d  were 
found,  the  process  terminated,  and  a  and  b  were  determined. 

For  a  charging  profile  at  any  age  of  life,  it  can  be  represented  as 
(U,  t ),  according  to  Eq.  (9),  with  all  parameters  fit  by  the  nonlinear 
least  squares  method.  As  the  battery  ages,  the  voltage-time  data 
sets  (U,  t)  shall  be  updated,  and  only  five  parameters  (i.e.  k.  At,  c,  d, 
e)  in  Eq.  (9)  shall  be  considered  in  the  iteration. 


U'(SOC)  =  (Joe fSOC)  -  (u;  +  / ■  RJ, )  exp  * soc  ■  exp5i  +  i(r'p+  R'q ) 

(5) 

According  to  Eqs.  (4)  and  (5),  the  voltage  diversity  with  variable 
SOC  is  expressed  as 

U'(SOC)  -  U(SOC)  =  (Up+I  Rp)exp-%  soc  expS 

-  (u'p  +rRp)exp“7ui'S0C-exp75i  (6) 

+  l(R>p  +  R>Q)-I{Rp  +  RQ) 

It  is  important  to  note  that  in  Eq.  (6),  Uoc(SOC)  is  not  present, 
which  eliminates  its  exhausting  requirement  for  measurement  in 


3.  Experimental 

Six  commercial  LIBs  (No.  1-6)  made  for  electronic  products 
were  tested  for  cycle  aging  experiments.  When  the  nominal  ca¬ 
pacity  faded  to  80%  of  the  rated  capacity,  the  test  was  terminated 
and  charge/discharge  data  analyzed.  The  test  schemes  are  formu¬ 
lated  as  follows: 

3.1.  Test  scheme  1 

Three  new  2.2Ah  LS-18650  (LiCo02/LiNixMnyCoz02)  batteries 
(No.  1-3)  were  charged  using  a  typical  CC-CV  method  where 
CC  =  C/2  with  an  end-of-charge  voltage  EOCV  =  4.2  V  and 
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CV  =  4.2  V  with  an  end-of-charge  current  EOCC  =  C/20.  In  the 
discharging  tests,  the  end-of-discharge  voltage  EODV  =  3.0  V  is 
used.  Battery  No.  2  was  subjected  to  60  °C  test,  while  No.  1  and  3 
were  tested  at  25  °C.  Battery  No.  3  was  subjected  to  test  with  a 
discharge  current  DC  =  3C  without  rest  interval,  while  No.  1  and  2 
were  discharged  at  1  C  with  10  min  rest  interval  between 
charge  and  discharge  regimes.  The  detailed  test  conditions  for 
each  battery  are  specified  in  Table  1.  As  Li  plating  may  occur  during 
low- temperature  charging  [12)  and  lead  to  additional  battery 
degradation,  low-temperature  tests  were  not  considered  in  this 
work. 

3.2.  Test  scheme  2 

Battery  No.  4  (LS-18650),  No.  5  (EVE-402030)  and  No.  6  (TJ-18650) 
with  different  cathode  compositions  and  capacity  ratings  were 
tested  at  room  temperature  using  conditions  listed  in  Table  2.  All 
subjected  to  CC  =  1C  and  DC  =  1C,  and  a  rest  interval  of  10  min  was 
applied  between  charge  and  discharge  regimes. 

4.  Results  and  discussion 

4.1.  Model  extraction 

In  Test  scheme  1,  the  benchmark  charging  profile  at  100%  SOH 
was  produced  from  the  first  cycle  of  battery  No.  1,  and  its  function  is 
denoted  as  U(t) lsi-  For  battery  No.  4-6  in  Test  scheme  2,  the 
benchmark  charging  profile  functions  are  denoted  respectively  as 
U{t) ls2,  D(t)EVE,  and  U(t) tj.  Four  function  formula  were  extracted, 
following  the  procedures  described  in  Refs.  [13],  as  follows: 

U(t) LSI  =  —0.30  exp-18  77-t  +  0.42-t - 0.35  t2  +  0.15- 13  +3.71 

(10) 

l/(f)LS2  =  -0.35  exp~46  02  r  +  1.22  f- 2.20- 12  + 1.81  f3  +  3.70 

(11) 

U(t) ^  =  -O.lO  exp-113  00  f +  0.83- 1  - 1.35- 12  +  1.20- 13  +  3.78 

(12) 

U(t)T}  =  -0.55  exp-4 10  t 1.12-t  +  1.30-t2  -  0.22-t3  +4.10 

(13) 

Using  the  method  described  in  Section  2,  parameters  a  and  b  in 
Eq.  (9)  were  identified  by  fitting  the  benchmark  charging  profiles 
respectively  for  each  battery.  The  SOH  estimation  using  the  math¬ 
ematical  model,  denoted  as  U'(t') lsi,  U'(t' )ls2 ,  U'(t') eve,  and  U'(t ')tj 
respectively,  can  be  extracted  as  shown  below: 

U'(t') LSI  =  -0.30-exp_18-77‘k(t+Atl  +  0.42  k(t'  +  At) 

-  0.35  •  k2{t’  +  At)2  +  0.1 5  •  k3  (f'  +  At)3  +3.71 
+  0.52-exp-2'47'k(t,+At)  -c  expd  d  +e  (14) 


Table  1 

Partial  specifications  of  Test  scheme  1. 


No. 

Types 

Cathode 

Capacity 

CC 

DC 

Int.  time 

Temp. 

1 

LS-18650 

LiCo02/LiNixMnyCoz02 

2200  mAh 

0.5  C 

1  C 

10  min 

25  °C 

2 

LS-18650 

LiCo02/LiNixMnyCoz02 

2200  mAh 

0.5  C 

1  C 

10  min 

60  °C 

3 

LS-18650 

LiCo02/LiNixMnyCoz02 

2200  mAh 

0.5  C 

3  C 

0 

25  °C 

Table  2 

Partial  specifications  of  Test  scheme  2. 


No.  Types 

Cathode 

Capacity 

CC 

EOCV 

EOCC 

EODV 

4  LS-18650 

LiCo02/LiNixMnyCoz02 

2200  mAh 

1C 

4.2  V 

0.05  C 

3.0  V 

5  EVE-402030 

LiCo02 

190  mAh 

1C 

4.2  V 

0.02  C 

3.0  V 

6  TJ-18650 

LiNixMnyCoz02 

1000  mAh 

1C 

4.3  V 

0.05  C 

2.7  V 

U'(t')  ls2  =  -0.35-exp_4602  k^+At)  +  1.22-k(t'  +  At) 

-  2.20-/<2(t'  +  At)2  +  1.81  -k3(t'  +  At)3  +  3.70 

+  0.52exp~1  •31k(t'+Af)  -c  expd  t  +e  (15) 

U'(t')  eve  =  -0.10-exp_11300  k(t,+At)  +  0.83  fc(t'  +  At) 

-  1.35-k2(t'  +  At)2  +  1.20-/<3(t'  +  At)3  +3.78 

+  0.57-exp~a21'k(t,+At)  -c  expd  d  +e  (16) 

l/'(t')TJ  =  -0.55-exp~4-1°'k(t,+Atl  -  1.12-fc(t' +  At) 

+  1.30-k2(t'  +  At)2  -  0.22  -l<3(t'  +  At)3  +4.10 
+  0.91  •exp-001'k(l,+A<l  —  cexpd  1  +e  (17) 

4.2.  Validation  analysis 

The  SOH  estimation  model  expressed  by  Eq.  (14)  is  used  for 
batteries  No.  1-3  to  estimate  the  SOH  through  cycle  aging,  respec¬ 
tively.  The  error  in  estimation  (SOHest.)  against  experimental  data 
(SOHeXp.)  is  illustrated  in  Fig.  3.  The  maximum  error  (discrepancy 
between  estimation  and  experiment,  ASOH  =  SOHest.  -  SOHexp.)  is 
mostly  below  2%.  For  batteries  No.  4-6,  with  each  estimation  model 
shown  in  Eqs.  15-17  respectively,  the  results  are  shown  in  Fig.  4,  with 
the  maximum  error  below  3%. 

There  are  additional  aspects  that  are  worth  noting  here.  The 
noticeable  fluctuations  in  error  (some  >2%)  around  800-950  cycles 
in  Fig.  3(a)  were  likely  caused  by  excessive  noise  from  the  experi¬ 
mental  data,  not  by  the  model  inaccuracy.  The  origins  of  these 
fluctuations  are  not  clear.  However,  it  is  worth  noting  that  after 
such  a  period,  the  trend  line  of  the  error  function  settled  around  1% 
level,  breaking  away  from  the  original  trend  of  increase.  Additional 
analysis  is  in  progress  to  investigate  if  this  transition  in  the  trend 
line  has  any  physical  implication  due  to  causality  from  the  test 
conditions  or  battery’s  properties.  Fig.  3(b)  shows  the  model  fol¬ 
lows  experimental  data  well,  including  a  rapid  descend  in  SOH  in 
the  vicinity  of  500  cycles.  Such  an  illustration  is  a  favorable  sign  to 
show  the  viability  of  the  model  prediction.  However,  one  should  be 
cautioned  to  use  this  sensitivity  in  predictability  mindfully  to  avoid 
being  deceived  by  artifacts.  In  Fig.  3(c),  once  the  SOH  falls  below 
80%,  the  capacity  fade  fluctuates  intensively;  so  is  error  function.  In 
this  region,  ASOH  is  close  to  4%  and  seems  to  increase  progressively. 
The  reason  is  not  clear  and  additional  investigation  is  under  way. 

The  examples  in  Fig.  4  illustrate  the  model  is  capable  of  esti¬ 
mating  SOH  in  different  cell  designs  and  chemistries  with  a 
reasonable  level  of  error  spreading  around  2—3%.  In  these  cases,  the 
charging  rate  has  been  increased  from  C/2  used  in  the  previous 
samples  to  1C.  Although  it  is  marginal  to  imply  if  a  higher  charging 
rate  such  as  1 C  would  lead  to  such  an  increase  in  error  from  2%  to 
3%,  it  is  nonetheless  worth  noting  that  higher  rate  charging  should 
increase  polarization;  thus  noisier  conditions  in  testing,  and  as  the 
rate  increases,  it  may  result  in  rate-dependent  causality  in  the 
disparity  of  state  functions  between  the  two  electrodes  in  the  cell. 
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Fig.  3.  The  SOH  estimation  against  experimental  data  and  the  error  (A SOH  =  SOHest  -  SOHexpJ  in  each  cycle  estimated  by  using  Eq.  (14).  Plots  (a)— (c)  are  for  batteries  No.  1—3 
respectively. 


Fig.  4.  The  SOH  estimation  against  experimental  data  and  the  error  (A SOH  =  SOHest  -  SOHexpJ  in  each  cycle  was  calculated  by  Eqs.  15-17,  respectively.  Plots  (a)-(c)  are  for  batteries 
No.  4-6,  respectively. 
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Using  graphite  as  the  negative  electrode  as  an  example,  the  high 
rate  charging  could  create  rate  limitation  situations,  or  worse 
lithium  deposition,  on  the  negative  electrode  [14]  and  lead  to  a 
more  severe  degradation  in  cell  performance.  Such  an  influence  on 
the  accuracy  of  the  model’s  predictability  cannot  be  ignored.  We 
thus  suggest  that  a  proper  charging  current  (e.g.  0.5C)  should  be 
selected  and  preferred  for  SOH  estimation,  as  the  lower  the  rate  the 
less  the  error  and  the  better  reliability  should  be  received.  Although 
we  have  not  discussed  the  influences  of  aging  and  degradation 
mechanisms  on  the  accuracy  of  the  model  prediction,  it  is  a  subject 
will  be  elaborated  in  our  future  work. 

5.  Conclusion 

An  SOH  estimation  method  is  illustrated  here,  which  is  based  on 
a  simple  equivalent  circuit  modeling  approach  to  parameterize  a 
single-variable,  time-based  SOH  inference  model  using  CC  charging 
profiles  at  various  stages  of  life  for  the  SOH  estimation.  The  initial 
CC  charging  profile  when  a  battery  is  commenced  for  service  is 
benchmarked  as  100%  SOH.  Subsequently,  the  SOH  of  the  battery 
can  be  estimated  by  this  inference  model  as  the  battery  ages.  The 
model  parameters  can  be  derived  from  a  nonlinear  least  squares 
algorithm.  The  high  accuracy  of  this  SOH  estimation  approach  was 
validated  by  experimental  results  using  cells  under  different 
charging  protocols  and  using  cells  with  different  sizes  and  chem¬ 
istries.  The  simplicity  of  this  method  makes  implementation  easy 
for  on-board  applications,  enjoying  the  benefits  of  less  demand  on 
prior  learning  and  universality  in  application  to  different  battery 
designs,  chemistries,  and  operating  conditions.  Especially,  this 


approach  gives  accurate  SOH  estimates  irrelevant  to  duty  sched¬ 
ules,  making  this  approach  attractive  for  practical  applications. 
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